Surface Josephson plasma waves in layered superconductors 
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We predict the existence of surface waves in layered superconductors in the THz frequency range, 
below the Josephson plasma frequency ujj. This wave propagates along the vacuum-superconductor 
interface and dampens in both transverse directions out of the surface (i.e., towards the super- 
conductor and towards the vacuum). This is the first prediction of propagating surface waves in 
any superconductor. These predicted surface Josephson plasma waves are important for difi'erent 
phenomena, including the complete suppression of the specular reflection from a sample (Wood's 
anomalies) and a huge enhancement of the wave absorption (which can be used as a THz detector). 

PACS numbers: 74.72.Hs 



High temperature Bi2Sr2CaCu208+5 superconductors 
have a layered structure of superconducting Cu02 planes 
with Josephson coupling between them [1]. This struc- 
ture favors the propagation of electromagnetic waves, 
called Josephson plasma waves ( JPW) [2-4] , through the 
layers. These waves attract considerable interest because 
of their Terahertz frequency range (from 300 GHz to 30 
THz corresponding to 1000-90 /im (wavelength)), which 
is still hardly reachable for both electronic and optical de- 
vices [6]. During the last decade there have been many 
attempts to push THz science and technology forward 
[6] because of many important applications in physics, 
astronomy, chemistry, biology and medicine, including 
THz imaging, spectroscopy, tomography, medical diag- 
nosis, health monitoring, environmental control, as well 
as chemical and biological identification. 

The unusual optical properties of layered superconduc- 
tors, including reflectivity and transmissivity, caused by 
the JPW excitation were studied in, e.g., Refs. [2, 7-10]. 
In particular, Ref. [8] demonstrated that the spectrum 
of JPW consists of two branches due to a peculiar ef- 
fect [11] of dynamical breaking of charge neutrality and, 
therefore, the transmissivity exhibits a sharp peak at fre- 
quencies just above the Josephson plasma frequency toj. 

All previous works on this problem have focused on 
running waves in the frequency range above the gap of 
the JPW spectrum, i.e., above the Josephson plasma fre- 
quency uj > ujj. A similar gapped spectrum also appears 
in solid state plasma [12]. In such situations the pres- 
ence of the sample boundary can produce a new branch 
of the wave spectrum inside the gap, i.e., a surface plas- 
mon [12, 13]. In general, surface waves play a very im- 
portant role in many fundamental resonant phenomena, 
such as the "Wood anomalies" in the reflectivity [13, 14] 
and transmissivity [15] of periodically corrugated metal 
and semiconducting samples, and are employed in many 
devices. Therefore, it is important to know if surface 
waves can exist in layered superconductors. 

Here, we prove that surface Josephson plasma waves 



can propagate along the surface between the supercon- 
ductor and the vacuum in a wide frequency range below 
LUj. We derive the dispersion relation for these waves and 
propose ways to excite these. We show that these sur- 
face waves play an important role in the absorption and 
reflection of electromagnetic waves, including their reso- 
nance dependence on the incident angle 9. The studied 
resonant absorbability could be experimentally observed 
by measurement of the surface impedance of a sample or 
dc resistivity. The predicted phenomena are potentially 
useful for designing THz detectors. 

Surface Josephson plasma waves. — We consider a 
semi-infinite layered superconductor in the geometry 
shown in the inset in Fig. 1. Using a standard approx- 
imation, the spatial variations inside the very thin su- 
perconducting layers are neglected in the direction per- 
pendicular to the layers. The crystallographic ab-plane 
coincides with the xy-plane and the c-axis is along the 
z-axis. Superconducting layers are numbered by the in- 
teger I >0. 

We consider a surface p-wave with the electric, E = 
{Ex,0, Ez}, and magnetic, H = {0, H, 0}, fields damped 
away from the interface z = 0, 

H, Ex, Ez oc exp(— iwt -I- iqx — klD) (1) 

inside a layered superconductor, z < 0, and 

^vac^ £;™=, El"" OC expi^iujt + iqx - Ke.eZ), (2) 

in vacuum, for z > and q> lo/c. Here D is the spatial 
period of the layered structure. 

The Maxwell equations for waves (2), qH^"'^ — 

-(w/c)^;r^ fcvacif™" = -(^w/c)^;™^ iqEi^^ + 

kvacE^"'^ — —{iuj/c)H'"^'^, provide both the usual disper- 
sion relation, fcvac — \/ (p' — I , for the wave in vac- 
uum and the ratio of amplitudes for tangential electric 
and magnetic fields at the interface z = -1-0 (i.e., right 
above the sample surface) 
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Inside the layered superconductor, where ^ < 0, the 
gauge invariant phase difference is described by a set of 
coupled sine-Gordon equations. For Josephson plasma 
waves, these non-linear equations can be linearized and 
rewritten in terms of the magnetic fields i?' between the 
Ith and {I + l)th layers, 
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Here Xab is the London penetration depth in the z- 
direction, the operator df is defined as df fi = + 
fi-i — 2fi, uij — y^STTeDJc/he is the Josephson plasma 
frequency determined by the maximum Josephson cur- 
rent Jc and dielectric constant e, and tOr is the relaxation 
frequency. The effect of breaking charge neutrality [11], 
which is crucial for our analysis, is taken into account in 
Eq. (4). The constant a characterizing this effect was 
estimated, e.g., in Ref. [9], a ^ 0.05-0.1 for Bi-2212 or 
Tl-2212 crystals. 

Substituting the wave (1) in Eq. (4), we obtain the 
implicit equation for the damped-wave transverse wave 
vector k{q,Lu) 
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layer. We ignore the displacement current in Eq. (6) be- 
cause it is proportional to a small parameter (Aabw/c)^. 
Besides, we neglected the contribution of the gauge- 
invariant scalar potential into the x-component of the 
electric field since the Debye length is much shorter than 
the wave length q^^ . In order to eliminate H", we use the 
relation i?™'^ - H° ^ i7™'=[l - eyip{-kD)], that follows 
from Eq. (1) with I = and / = 1. Now, using Eq. (6), 
we obtain the ratio between electric and magnetic fields 
at z = —0 (i.e., right below the sample surface) 
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(7) 



Using the continuity conditions for the tangential com- 
ponents of electric field at the surface and Eqs. (3), 
(5), and (7), we obtain the dispersion relation for two 
branches of the surface wave corresponding to two solu- 
tions of Eq. (5). For {\—uj/u)j) » ^ajeDjXah ~ 
5 ■ 10^"*, this spectrum can be written as 



x{l + 2Tn [l-(l + r,V)'^']}''^' 



(8) 



for two branches 
sionless variables: 



■'=F" . Here we introduce the dimen- 

K = cq/ujj, fl = w/ujj, (3 — 



2Xl^ujj/cD, and Fn = (1 - ri^)/4a. The value of the pa- 
rameter (3 for Bi-2212 is about 1.4. The spectra Eqs. (8) 
are shown in Fig. 1. Both branches merge in a narrow fre- 
quency region below wj, i.e., (1—0) ~ {a/ey^^{D/Xab)- 



for LUr 0. A similar spectrum, but for running Joseph- 
son plasma waves (with imaginary k[q,uj)), was earlier 
obtained in Ref. [9]. At w < wj, solving Eq. (5) with 
respect to sinh^[fc(g, a;)£'/2] results in two branches of 
positive transverse spatial decrement, fc±(g,aj) > 0. It 
is important to stress that Eq. (5) is not a spectrum of 
the studied longitudinal surface waves. Indeed, Eq. (5) 
has two free parameters, u) and q. We need to obtain 
the q{oj) dispersion relation for our surface waves. This 
dispersion relation q{uj) can be obtained by joining fields 
(1) in a superconductor and (2) in vacuum at the sample 
surface via the boundary conditions. Thus, in order to 
find the spectrum of the surface Josephson plasma waves, 
we can derive the ratio Ex/ H in the superconductor at 
the sample surface and equate this ratio to Eq. (3). 

The difference between the magnetic field H^^'^ in vac- 
uum and the value iJ" between the 0th and the 1st super- 
conducting layers is described by the London equation. 
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Here Axo and Exo are the x-components of the vector 
potential and electric field in the first superconducting 




FIG. 1: (Color online) Spectra of the surface Josephson 
plasma waves for the parameters a = 0.1, and P — 1.4, 
standard for Bi-2212. The inset shows the geometry used. 



The surface mode K_(ri) attenuates into the vacuum 
at very large distances, of about cXat/i-^D, and seems to 
be difficult to observe experimentally. Another mode. 
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K+(f7), dampens on scales ^ cjuj and is of significant 
interest. As we show below, it plays an important role in 
the transmissivity and reflectivity properties of layered 
superconductors at frequencies uj <loj. 

Excitation of the surface waves; resonant electromag- 
netic absorption. — One of the ways to excite surface 
waves is via externally applied electromagnetic waves on 
a sample having spatially modulated parameters. Thus, 
we consider a weak modulation of the maximum current 
density Jc via, say, creation of pancake vortices by the 
out-of-plane magnetic field. This can result in the mod- 
ulation of wj. For simplicity, we assume that 
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M « 1, (9) 



where a is a spatial period. 

An electromagnetic wave with uj < ujj incident at 
an angle 6 with respect to the sample surface gen- 
erates modes having longitudinal wave vectors q„i = 
wsin^/c + 27rm/a, with integer m. Almost all of these 
modes for m ^ are weak, because jj, 1. However, 
one of these modes (e.g., for m = 1) can be excited with 
large amplitude at resonance, i.e., when the wave vec- 
tor qi = {u)/c)sm6 + 2Tr/a is close to the wave vector 
(7_i- = LOjK-^-{Cl)/c of the surface wave (8). This corre- 
sponds to the incident angle 9 close to the resonance an- 
gle ^0 defined by 

2'KC 

f2sin6'o-h = K+(f2). (10) 

OLOJ 

Because of this resonance, the amplitude of the wave with 
q = qi can be of the same order, or even higher than, 
the amplitude i?™*^ of the incident wave. In resonance, 
the mode with q = qi is actually the surface Josephson 
plasma wave discussed above. 

Let us discuss the mechanism of excitation of surface 
waves. The incident -ff;™'^ (red arrow in Fig. 2) and 
specularly-reflected H^^^ (blue dashed arrow in Fig. 2) 
electromagnetic waves 
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(11) 



generate the wave (in green. Fig. 2) damped inside the 

superconductor, 

Hoex.^U[ — ijw-/t_( — - — ,uj\Dl\, (12) 

having a lower «;_(0), compared to k+{Q.). Indeed, as 
was shown in Ref. [9], the longer wave is mainly gen- 
erated at (w — wj)/wj » oijeD jXab ~ 5 • 10~^. 
Moreover, at low frequencies, this longer wave exhibits a 
correct limiting behavior, producing magnetic fields de- 
scribed by the well-known London equation. The waves 
Hy^, Hl^ « Ry^, and w 2^;^^'^ represent the solu- 
tion of the problem in the zero approximation with re- 
spect to \i. 



To first order approximation, the solution of Eq. (4) 
with u)j{x) inside the superconductor consists of both 
the forced component. 



iJforced exp < iq\x — itot — K- 



wsmt 



Lo]Dl} , (13) 



which attenuates with the same decrement as Hq, and 
free oscillations 

i?surf exp {iqix — icot — k+ {qi ,ui)Dl} . (14) 

The amplitude ifforced is defined by Eq. (4) itself, whereas 
the amplitudes, ffsurf (in dotted magenta) and corre- 
sponding vacuum mode 

H^^ exp {iqix-iujt- (ql-^] z\ , (15) 



are determined by the boundary conditions for harmonics 
q = qi oi H and aX z = 0. Solving the correspond- 
ing set of nonhomogeneous linear equations for ilsurf and 
-^surf ^6 obtain the amplitude of the resonance wave: 



TTvac _ rrvac 
surf in 



with the denominator 
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(16) 



■R = X + iY = 



2k+(0) cos6»o(6'-6io) 
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(17) 

The resonance is characterized by the Lorentz form 
1/TZ = l/[X{9) + iY]. The excitation of the surface wave 
results in the resonant peak in the electromagnetic ab- 
sorption. 



Absorption(6') oc a±E'^ + anEl cx 



X2(6i)-Fy2 



, (18) 



shown in Fig. 2. The resonance in the absorption can 

be observed by mcasiiring the dependence of the surface 
impedance on the angle 9. Alternatively, the peak in ab- 
sorption produces a temperature increase, resulting in a 
sharp increase of the DC resistance or even the transition 
of the sample to the normal state a.t 9 = 9o. 

Naturally, the absorption peak is accompanied by the 
resonant decrease of the amplitude H^p^ of the specu- 
larly refiected wave. Even though this effect is of second- 
order with respect to /z, it can result in a significant (up 
to complete) suppression of -ff™*^ at 6* = 6*0 due to the 
resonance denominator TZ{9). In the optical range, this 
phenomenon is known as Wood's anomalies [14] of reflec- 
tivity and has been used for several diffraction devices. 

THz detectors. — These effects could be potentially 
useful for the design of THz detectors, an important cur- 
rent goal of many labs worldwide. The simplest design 
could be a spatially modulated Bi2212 sample fixed on 
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incident wave specular wave 




FIG. 2: (Color online) (a) Schematic diagram showing the 
mechanism of excitation of the surface waves along the 
superconductor-vacuum interface. To zero-order approxima- 
tion, with respect to the amplitude /i of the spatial modula- 
tions in a superconducting sample, an incident wave (shown 
as a solid red arrow) reflects as a specular wave (the straight 
dashed blue arrow), producing a damped wave (green wave) 
inside superconductors. To first order approximation, the 
very intense surface wave (dotted magenta wave) can be 
excited at a certain resonant angle between the incident 
wave and sample surface, (b) Absorption obtained using 
Eq. (18) for different effective dampings 5 = /3a;rQ!fi*(fc+ — 
^2)1/2 [2ljj(^i _ J72)3fc+ cos6lo]"^ 

a precisely rotated holder and attached by contacts to 
measure its resistance. Spatial modulations in the sam- 
ple could be fabricated by either using ion irradiation of 
the sample covered by periodically modulated mask [16], 
or even mechanically [17]. When rotating the sample, 
the incident THz radiation can produce a surface wave 
at certain angles. This results in a strong enhancement 
of absorption associated with increasing of temperature 
in the sample and, thus, its resistance. The relative po- 
sitions of the resonance peaks (the set of angles) allows 
to calculate the angle and the frequency of the incident 
THz radiation, while the relative heights of the resistance 
peaks can be used to estimate the intensity of the inci- 
dent radiation. 

Conclusions. — We have derived the surface Joseph- 
son plasma waves in layered superconductors and ob- 



tained their dispersion relation. The absorption of the 
incident electromagnetic wave can strongly increase at 
certain incident angles due to the resonant generation of 
the predicted surface waves. This is the first prediction 
of propagating surface waves in any superconductor. We 
propose a way to experimentally observe these surface 
waves. 
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